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We study the effects of magnetic and electric fields on the g factors of spins confined in a two- 
electron InAs nanowire double quantum dot. Spin sensitive measurements are performed by moni- 
toring the leakage current in the Pauli blockade regime. Rotations of single spins are driven using 
electric-dipole spin resonance. The g factors are extracted from the spin resonance condition as a 
function of the magnetic field direction, allowing determination of the full g-tensor. Electric and 
magnetic field tuning can be used to maximize the ^-factor difference and in some cases altogether 
quench the EDSR response, allowing selective single spin control. 

PACS numbers: 85.35. Gv, 71.70.Ej, 73.21.La 



The desire for a fully controllable and scalable quan- 
tum computer places several stringent conditions on its 
constituent qubits [Tjl. Recent implementations of spin- 
qubits in GaAs heterostructures have made impressive 
advances, resulting in a well characterized system in 
which spin initialization, control, and readout have all 
been achieved [IHl]. For GaAs-based qubits, rapid and 
selective control of single spins remains a challenge. Sin- 
gle spin rotations have been performed using traditional 
electron spin resonance (ESR) as well as electric-dipole 
spin resonance (EDSR), but the highest achieved manip- 
ulation rate in GaAs is still 2 orders of magnitude slower 
than the exchange gate [2l |3l El |6] . Encoding the qubit in 
two- or three-electron spin states can eliminate the need 
for single spin control, but requires local magnetic field 
gradients or exquisite control of the exchange interaction 

[ail]- 

"Spin-orbit qubits" based on InAs nanowires have re- 
cently been proposed as an alternative to the GaAs lat- 
eral quantum dot system [9^^. In the spin-orbit qubit, 
the qubits are dressed states of spin and orbital degrees 
of freedom, due to the strong spin-orbit coupling of InAs. 
Electrical control of the qubit 's orbital component allows 
Rabi frequencies on the order of 60 MHz to be achieved 
m Ho] . While strong spin-orbit coupling enables fast spin 
rotations, it has the potential to introduce several com- 
plications. In particular, the electronic g factors, spin 
relaxation time, and EDSR rotation rates are expected 
to vary with magnetic field direction [11 . Moreover, the 
spin-orbit interaction couples the spin to the orbital com- 
ponent of the wavefunction, leading to a g-tensor that is 
sensitive to the gate-tunable confinement potential. 

In this Letter, we explore the effects of strong spin- 
orbit coupling on spins confined to an InAs nanowire dou- 
ble quantum dot. We demonstrate magnetic and electric 
field control of the g factors for each quantum dot. The 
EDSR spin manipulation rate is a sensitive function of 
magnetic field direction, allowing the Rabi frequency to 
be maximized. For specific magnetic field directions the 
EDSR response can be dramatically reduced. Our results 
show that electric and magnetic field tuning of the g fac- 



tor can be used to optimize the EDSR driving rate and 
allow for selective single spin control. 

Our device, shown in Fig. 1(a), is fabricated on a 
high resistivity, oxidized silicon substrate. Using elec- 
tron beam lithography, we define an array of gate pat- 
terns consisting of a series of 20 nm thick Ti/Au gates, 
spaced at a 60 nm pitch. Two large side-gates allow the 
transparency of the nanowire leads to be tuned. A 20 nm 
layer of SiN^^ is then deposited as a gate dielectric using 
plasma enhanced chemical vapor deposition [10 . Single 
crystal InAs nanowires with a zinc blende structure are 
grown using a gold catalyzed vapor- liquid- solid process 
in a metal-organic chemical vapor deposition reactor [T^ • 
Nanowires are removed from the growth substrate with 
ultrasonication in ethanol and then dispersed on the gate 
array. The final fabrication step involves defining low 
resistance ohmic contacts to nanowires with diameters 
of ~50 nm that have fallen across a gate pattern. The 
sample was measured in a dilution refrigerator equipped 
with a vector magnet system and high frequency coaxial 
wiring. 

We first determine the charge stability diagram by 
measuring the current at finite applied bias, as shown in 
Fig. 1(b). The absence of finite bias triangles in the lower 
left hand corner of the plot indicates the double quantum 
dot has been completely emptied of free electrons. In this 
region (Nl^Nr) = (0,0), where Nl{Nii) is the number 
of electrons in the left (right) dot. Measurements at high 
bias and with relatively transparent tunnel barries have 
confirmed our identification of the (0,0) charge state. We 
focus on the two-electron regime, where Pauli blockade 
results in current rectification at the (1,1)^(2,0) charge 
transition, as first observed in vertical double quantum 
dots by Ono et al [14]. At positive bias, charge transport 
at the (l,l)<->(2,0) transition occurs freely through a cy- 
cle of (1,0)^(2,0)^(1,1)^(1,0), where steps with double 
occupancy are limited to spin singlet states due to the ~9 
meV exchange splitting of the (2,0) state. Under nega- 
tive bias, the cycle is reversed and the system loads the 
second electron into the (1,1) charge state, where the sin- 
glet and triplet states are nearly degenerate. If a (1,1) 
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FIG. 1: (Color online) (a) Scanning electron microscope im- 
age illustrating the device geometry. Four gates are sufficient 
to fully define the double quantum dot with this device, re- 
sulting in a double well confinement potential. The double 
dot's coupling to the right and left leads are controlled with 
voltages Vrw and Viw respectively, while Vrp and Vip control 
the interdot coupling and occupation, (b) Current measured 
through the device as a function of gate voltages Vrp and Vip , 
with source drain bias Vsd — -4 mV. Dashed lines are super- 
imposed to illustrate the double dot charge stability diagram. 
Left inset: (1,1) ^ (2,0) transition at Vsd = +4 mV. The 
apparent excited state is a resonance associated with the one- 
dimensional leads that contact the double quantum dot [12]. 
Right inset: (1,1) ^ (2,0) transition at Vsd — -4 mV. The 
current in the insets has been multiplied by a factor of two. 
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FIG. 2: (Color online) Leakage current measured as a func- 
tion of magnetic field and applied frequency for (a) B 1 1 z and 
(b) B II [x -\- z). The zero field peak, shown in (c), is due to 
hyperfine field mixing of the spin states, with a fit to the data 
yielding Bn — 3.3 mT. The spin resonance condition is sat- 
isfied when hf = ^/iB|B|, resulting in a peak in the leakage 
current. Cuts through the data in (a) and (b) are shown in 
(d), demonstrating the dependence of the g factors on mag- 
netic field direction. The upper trace is offset by 0.3 pA for 
clarity. 



triplet state is loaded, the charge transition to (2,0) will 
be blocked due to Pauli exclusion. The left and right 
insets to Fig. 1(b) show the current through the double 
dot at positive and negative bias, exhibiting the voltage 
bias dependence characteristic of Pauli blockade. 

In the absence of spin-flips or cotunneling, the leakage 
current in the Pauli blockade regime will be zero. Pro- 
cesses which drive rotations from the (1,1) triplet states 
to the (1,1) singlet state will lift Pauli blockade, resulting 
in a measurable leakage current. Pauli blockade enables 
spin-dependent readout and has been used to detect the 
mixing of spin states due to hyperfine fields, ESR, and 
EDSR in GaAs and InAs quantum dots P l6l [TOl [T5] . 

We measure the leakage current in the Pauli blockade 
regime as a function of magnetic field B and the fre- 
quency / of the ac excitation applied to a depletion gate 
(see Fig. 2). A zero field peak is evident, corresponding 
to rapid singlet-triplet mixing due to the hyperfine fields 
[151 [16]. We extract a hyperfine field B^ = 3.3 mT from 
these data using the theory of Jouravlev et al. [17 . Pre- 
vious reports indicate B^ = 0.66 and 1.5 mT in InAs 
quantum dots, with a value that depends on the number 
of nuclei the electron spin has appreciable overlap with 
[inillH]. Away from zero field, the leakage current is non- 
zero only when the EDSR resonance condition is satisfied. 
The data in Fig. 2(a) are measured for B || i, while the 
data in Fig. 2(b) are measured for B || {x -\- z). Vertical 
cuts through the data, shown in Fig. 2(d), clearly indi- 
cate that the g factor of each electron is sensitive to the 



magnetic field direction, with well-resolved EDSR peaks 
visible for each spin in the lower trace. 

We quantitatively measure the directional dependence 
of the g factors by fixing the magnetic field magnitude at 
|B| = 80 mT and rotating the field direction. Data are 
shown in Figs. 3(a)-(c) for three different field rotations. 
In each case the position and intensity of the EDSR reso- 
nances evolve with the field direction. For each magnetic 
field direction we fit each resonance peak to a Gaussian 
and extract the g factor. Fitting is only performed for 
regions of the data where the resonance is clearly visible. 
After the g factors have been extracted, we perform a 
simultaneous curve fit of the data sets in Fig. 3(a-c) to 
a general model of a 3D anisotropic g factor: 
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Here , ^2 , and ^3 are the values of the g factors for the 
three principal axes, and 5i, B2 and B^, are the magnetic 
field components along these axes [I9l [20] . Our fitting 
also takes into account field offsets of 2-3 mT due to 
trapped flux in the superconducting magnet coils. 

The best fit values for the g tensor are listed in Ta- 
ble 1. The principal axes' g factors vary from approxi- 
mately 7 to 9, suppressed from the bulk value of |^| = 
14.7 through a combination of spin-orbit coupling and 
confinement [2TH23] . In each case, the g tensor contains 
one large axis, and two approximately equal and smaller 
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FIG. 3: (Color online) (a-c) Leakage current measured as a function of magnetic field direction and frequency for the 
rotations depicted in the insets below. The y-axes have been rescaled from frequency to units of g, taking a fixed magnetic 
field amplitude of 80 mT; all angles are radians. In this coordinate system, the unit vector describing the nanowire axis is 
approximately (^,0) = (1.2,3.1). A fit to Eqn. 1 has been superimposed on the data, (d-e) Pumped current extracted from 
the upper resonances in (a-c), along with a fit to Eq. (3). 



axes, appearing to mimic the cylindrical symmetry of the 
nanowire itself. However, we find that the large axis is 
oriented in different directions for each electron, and in 
general is not aligned with the nanowire axis. These re- 
sults suggest that the orientation of the g tensor is largely 
determined by the anisotropy of the confinement poten- 
tial, not the underlying one-dimensional symmetry of the 
nanowire. 

We performed a second set of measurements with the 
device tuned to a very different gate voltage configura- 
tion in order to measure the sensitivity of the g tensor to 
changes in the confinement potential. The data in Fig. 3 
were acquired with the double quantum dot tuned such 
that the tunnel rates to the left and the right leads were 
balanced. To explore a different confinement potential, 
we unbalanced the tunnel rates by varying the voltages on 
the local electrodes, thereby tuning the double quantum 
dot into a regime where transport shows a large, asym- 
metric cotunneling peak [see Fig. 4(a)]. Figure 4(b)-4(d) 
map out the g factor as a function of field orientation for 
the 'unbalanced' case. The g tensor has been dramati- 
cally altered, with extracted values for the principal axes 
directions and g factors summarized in Table I. For this 
gate voltage configuration, the upper resonance is highly 
anisotropic and shifted to higher frequencies. The lower 
resonance is broadened significantly, which we attribute 
to enhanced dynamic nuclear polarization [24 . 

We now focus on the field dependence of the EDSR am- 
plitude, which is expected to depend on both the mag- 
nitude and direction of the magnetic field [5]. Figure 
3(d)-3(f) show the current extracted from the upper res- 
onance of Fig. 3(a)-3(c) respectively. We find that the 
on resonance current is strongly modulated by field ori- 
entation. Golovach et al. |5j have shown that the ef- 



fective magnetic field for spin-orbit driven EDSR in a 
two-dimensional system is given by 

B,^(t) = 2B X (l^osincj^c^) (2) 

where ^Iq is a dimensionless "spin-orbit vector" deter- 
mined by the applied electric field and the local spin- 
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FIG. 4: (Color online) (a) Measured current through the 
double dot as a function of Vip and Vrp for Vsd — -4 mV in the 
case of unbalanced tuning, (b)-(d) Evolution of the resonance 
condition for magnetic field rotations about the x, y and z 
axes respectively. In (b)-(d) the y axes have been rescaled 
into units of g^ taking a fixed magnetic field amplitude of 80 
mT. 
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TABLE I: C-tensor Parameters. Values are given for both 
the balanced case of Fig. 3, and the 'unbalanced' case of Fig. 
4. a, /3 and 7 are Euler angles (in radians) which define the 
transformation from the measurement coordinate system to 
the g tensor coordinate system, using the ZX' Z" convention. 



Dot 


1^1 1 


1^2 1 


1^3 1 


a 




7 


Balanced ^ 
B 


9.1 
8.4 


7.8 
7.3 


7.5 
7.0 


1.9 
-0.81 


2.1 
1.0 


-0.25 
1.5 


Unbalanced 

B 


22 
8.8 


12 
7.6 


8.0 
7.4 


1.8 
-1.2 


1.9 
1.0 


-0.21 
0.73 



orbit parameters. In the absence of an applicable the- 
ory, we generalize the concept of the spin-orbit vector 
to three dimensions. While Eqn. 2 gives the Rabi fre- 
quency, fr = giiBBso/'^h for on resonance driving, the 
dependence of the pumped current is less clear. In gen- 
eral, the current will not be proportional to fr as the 
fluctuating nuclear field shifts the applied excitation on 
and off of resonance. We have studied the induced cur- 
rent using a simple rate equation model, similar to the 
approach used by Koppens et al. [25 . For Bso < ^at we 
find 
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where F^ is the interdot tunneling rate. Fits to Eqn. 3 
are displayed in Fig. 3(d)-3(f). Equation 3 reproduces 
the periodic modulation of the pumped current; however 
we are not able to simultaneously fit the three individual 
rotation sweeps. The discrepancy from the behavior pre- 
dicted by Eqs. 2 and 3 is not well understood but may 
involve corrections due to the magnetic field dependence 
of spin-orbit matrix elements [11 . 

Equation 3 may also be used to estimate the magnitude 
of the spin-orbit field. With a slightly different tuning, 
we have observed EDSR currents exceeding 1.5 pA at 
\B\ = 150 mT. For this tuning, we measure positive- 
bias (nonblockaded) currents of ~35 pA. Assuming the 
interdot tunnel coupling is rate limiting, this results in 
a lower bound estimate of F^ 220 MHz, leading to a 
spin-orbit field |Bso| = 1.0 mT. Finally, we may estimate 
the spin-orbit length Igo as 

where |£^| is the applied electric field, rrie is the effective 
mass, and A is the level spacing [5 . We estimate \E\ 
as 1 mV/60 nm, and find A = 9 meV from finite bias 
measurements. With these parameters, we find Igo ~ 170 
nm. Previous studies have reported values for the spin- 
orbit length in InAs nanowires of 100-200 nm [26) |27 ] . 

The large variation of the g factors with electric field 
tuning suggests that g tensor modulation might play a 
significant role in electrically driven spin manipulation 
[28 . Assuming the change in the g factor is purely linear 
in electric field, we find a maximum dg/dVg ~ 0.03/mV. 



Taking g = 8 and assuming an applied magnetic field of 
80 mT, this yields an effective ac magnetic field of 0.3 
mT, which is of the same order as the calculated spin- 
orbit field. Our results suggest that g tensor modulation 
may be significant for specific electric and magnetic field 
configurations. 

We have demonstrated magnetic and electric field tun- 
ing of the electron g factors in an InAs nanowire dou- 
ble quantum dot, which have important implications for 
spin-orbit qubits. The ability to selectively control sin- 
gle spins may be obtained by tuning magnetic and elec- 
tric fields to (1) maximize the difference between the g 
factors of neighboring spins and (2) maximize the differ- 
ence in the response of neighboring spins. A linear array 
of quantum dots in an InAs nanowire could be used to 
implement a Loss-DiVincenzo style spin register, where 
neighboring spins are tuned out of resonance with electric 
field tuning of the g tensor. The combination of spectral 
selectivity with the localized nature of the on-chip gen- 
erated microwave fields may allow the development of 
scalable quantum dot arrays in InAs nanowires. 
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